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In this paper I present a detailed QCD analysis of twist-3 effects in the Wandzura-Wilczek (WW) 
approximation in deeply virtual Compton scattering (DVCS) observables for various kinematical 
settings, representing the HERA , HERMES, CLAS and the planned EIC (electron-ion-collider) 
experiments. I find that the twist-3 effects in the WW approximation are almost always negligible 
at collider energies but can be large for low Q 2 and smaller Xbj in observables for the lower energy, 
fixed target experiments directly sensitive to the real part of DVCS amplitudes like the charge 
asymmetry (CA). Conclusions are then drawn about the reliability of extracting twist-2 generalized 
parton distributions (GPDs) from experimental data and a first, phenomenological, parameterization 
of the LO and NLO twist-2 GPD H , describing all the currently available DVCS data within the 
experimental errors is given. 
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I. INTRODUCTION 

Hard, exclusive processes and amongst them deeply 
virtual Compton scattering (DVCS) in particular 0, S 
0,0 SEISES El EI; have emerged in rece nt y ears as 
prime candidates to gain a three dimensional |18| image 
of parton correlations inside a nucleon p!4llalllill7lll8| . 
This information is gained by mapping out the key com- 
ponent containing this three dimensional information, 
namely generalized parton distributions (GPDs). 

GPDs have been studied extensively in recent years 
USSSSSSSSIlSinilll since these distribu- 
tions are not only the basic, non-perturbative ingredient 
in hard, exclusive processes such as DVCS or exclusive 
vector meson production, they are generalizations of the 
well known parton distribution functions (PDFs) from 
inclusive reactions. GPDs incorporate both a partonic 
and distributional amplitude behavior and hence contain 
more information about the hadronic degrees of freedom 
than PDFs. In fact, GPDs are true two-parton correla- 
tion functions, allowing access to highly non-trivial par- 
ton correlations inside hadrons [l9|. 

In order to perform a mapping of GPDs in their vari- 
ables, experimental data from a wide variety of hard, 
exclusive processes is needed. Furthermore, in order 
to be able to properly interpret the data, the processes 
should be well understood theoretically. Therefore, one 
should start out exploring the theoretically "simplest" 
process i.e. the one with the least theoretical uncer- 
tainty. In the class of hard, exclusive processes this is 
DVCS (e(fc) +p(P) -► e(k')+p(P') +j(q')). The rea- 
son for this is the simple structure of its factorization 
theorem S S The scattering amplitude is simply 
given by the convolution of a hard scattering coefficient 
computable to all orders in perturbation theory with one 
type of GPD carrying the non-perturbative information. 
Factorization theorems for other hard exclusive processes 
are usually double convolutions containing more than one 
non-perturbative function. 



As with all factorization theorems, the DVCS scat- 
tering amplitude is given in this simple convolution 
form up to terms which are suppressed in the large 
scale of the process. In this case the large scale 
is the transfered momentum i.e. the virtual photon 
momentum, Q, and the suppressed terms are of the 
order 0((m N /Q) n ), 0{{yf=t/Q) n , 0{(X QCD /Q) n ) with 
tojv the proton mass and t = (P — P') 2 the momentum 
transfer onto the outgoing proton. This means, however, 
that for smaller values of Q these uncontrolled terms, 
called higher twist corrections, could in principle be size- 
able and the convolution or lowest twist term need not 
be the leading one. Since we are interested in the extrac- 
tion of the non-perturbative information of the lowest 
twist GPD, in this case twist-2, we need to understand 
something about these higher twist corrections. 

It was shown by several groups [2(], |2l|, that 
the first suppressed term (twist-3) in the DVCS scatter- 
ing amplitude can, in leading order of the strong coupling 
constant (LO) and in the Wandzura-Wilczek (WW) ap- 
proximation, be simply expressed through a sum of terms 
which can be written as a convolution of a LO twist-3 
hard scattering coefficient with a twist-2 GPD. Unfor- 
tunately, it could not be shown that twist-3 in the WW 
approximation does factorize to all orders in perturbation 
theory (it seems to hold in NLO though 24]). However, 
we do have at least some control over the leading of the 
higher twist terms and can therefore use parameteriza- 
tions of twist-2 GPDs also for twist-3. 

Equipped with this knowledge one might think that 
this information is enough to obtain an unambiguous in- 
terpretation of DVCS data in terms of GPDs, but life is 
unfortunately even more complicated. Besides the dy- 
namical twist contributions to the amplitude, there are 
kinematical power corrections in the DVCS cross section 
i.e. contributions where a dynamical twist-2 amplitude 
is multiplied by a term oc (>/—i/Q, ttin/Q) which makes 
them of the same power in 1/Q as twist-3. These terms 
are of particular importance in the interference term be- 
tween DVCS and the QED Bethe-Heitler (BH) process 
which both contribute to the total DVCS cross section. 
Since the kinematical power corrections are nothing but 
dynamical twist-2 with a kinematical dressing they can 
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be handled by the same GPD parameterization as dy- 
namical twist-2. Thus we have the leading corrections to 
the DVCS process under control, at least in LO, and can 
now investigate several things: a) How big are physical 
observables at various values of Xbj , Q 2 and t for a certain 
center-of-mass energy s, especially those terms isolating 
various parts in the interference term directly propor- 
tional to real or imaginary parts of DVCS amplitudes? b) 
How big are the next-to-leading order (NLO) corrections 
in the strong coupling constant to these observables? c) 
How big are the corrections to these observables due to 
twist-3 effects? d) How reliably can the twist-2 GPDs be 
extracted from DVCS data? 

In this paper I will concentrate on c) and d). a) and 
b) have been extensively discussed using various GPD 
models in [25l I2H l27j . c) and d) have, in various forms, 
been discussed in |28ll2^ |. This was done, however, with- 
out taking evolution effects into account. We know that 
evolution effects are sizeable and can change the shape of 
the GPD in the ERBL region substantially [3(|, thereby 
strongly affecting the real part of the amplitude I3ll and 
the observables associated with the real part. In |25ll2(| 
it was demonstrated that the type of GPD models most 
commonly used cannot be brought into agreement with 
the available data once evolution effects were taken into 
account. Subsequently, in |32l |. a model was proposed 
which was able to describe all currently available DVCS 
data within a full NLO QCD analysis. 

I will extend the analysis of [321 ] , by not only using the 
input model of [32j for the twist-2 GPDs but also as an 
input for the twist-3 sector in the WW approximation. 
Sec. [H] contains the DVCS kinematics, structure of the 
cross section, equations for the twist-3 contributions and 
the definition of the relevant observables. In Sec. 11111 
I will recapitulate the model I use for the four relevant 
twist-2 GPDs H, H, E, E and their twist-3 counter parts. 
Sec. IIVI contains the twist-2 vs. twist-3 results for DVCS 
observables in various kinematical settings. In Sec. [3 I 
will present for the first time a phenomenological param- 
eterization of the twist-2 GPD H which can describe all 
currently available DVCS data within the experimental 
errors in both LO and NLO QCD. I will then summarize 
in Sec. |Vl| 



II. DVCS: KINEMATICS, CROSS SECTION 
AND DEFINITIONS 

The lepton level process, e ± (fc, n) N(P,S) — ► 
e ± (fc',K / ) N(P',S') j(q',e'), receives contributions from 
each of the graphs shown in Fig. This means that the 
cross section will contain a pure DVCS-, a pure BH- and 
an interference term. 

I choose to work in the target rest frame given in 
pEjj (see Fig. [5J , where the positive z-direction is chosen 
along the three-momentum of the incoming virtual pho- 
ton. The incoming and outgoing lepton three-momenta 
form the lepton scattering plane, while the final state 




FIG. 1: a) DVCS graph, b) BH with photon from final state 
lepton and c) with photon from initial state lepton. 




FIG. 2: The kinematics of the leptoproduction in the target 
rest frame. 



proton and outgoing real photon define the hadron scat- 
tering plane. In this reference frame the azimuthal angle 
of the scattered lepton is cf>i — 0, while the azimuthal 
angle between the lepton plane and the final state pro- 
ton momentum is <f>N = <j). When the hadron is trans- 
versely polarized (within this frame of reference) S±_ = 
(0, cos<P, sin<P, 0) and the angle between the polarization 
vector and the scattered hadron is given by (p — <l> — <j}^ . 
The four vectors are k = (E, E sm8i,0, E cosBi), q = 
(g°, 0, 0, — |g 3 |). Other vectors are P = (M, 0,0,0) and 
P' = (E', |P'|cos0sin6»Ar, sin0sin6>Ar, |P'| cos 6»jv). 
The longitudinal part of the polarization vector is Slp = 
(0,0, 0, A). The relevant Lorentz-invariant variables for 
DVCS are then: 



t-& 1 -V--t-*-*-<?->V-'-Ti- 

where P = (P + P')/2, q = (q + q')/2 and which are 
related to the experimentally accessible variables, C = 
xi,j = —q 2 /2P ■ q and Q 2 = —q 2 , used throughout this 
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paper, via 
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2-C(l-*) 2 ^ 
Note that t has a minimal value given by 



t- 



Q 



2 2(1 - Xbj ) (i - vtt?) + . 

4x6^(1 - Xbj) + e 2 



(1) 



(2) 



where e 2 = 4M 2 x 2 . j /Q 2 . Thus the theoretical limit f -> 
in an exclusive quantity is not attainable in any experi- 
mental set-up and one will have to rely on extrapolations. 
The corresponding differential cross section is given by 



The c n 's and s„'s are the Fourier coefficients of the 
cos(n0) and sin(n0) terms. These coefficients are given 
as combinations of the real and imaginary part of the un- 
polarized and the polarized proton spin-non-flip and spin- 
flip DVCS amplitudes H,H,S,E (for the c r 's or s x, s) or 
the squares of the afore mentioned DVCS amplitudes (for 



the c uvcb 's or s DVCb 's) . The exact from is given in |2 
and does not have to be repeated here. I will discuss 
the computation of the DVCS amplitudes and the neces- 
sary model assumptions in the next section. The precise 
form of the BH propagators Vx^{4>) which induces an ad- 
ditional (/(-dependence, besides the cos(n</>) and sin(n</>) 
terms, and which can mock cos(n0) and sin(n</f) depen- 
dences in certain observables, can also be found in |28| . 
Note that in order to avoid collinear singularities occur- 
ring through the incidence of the outgoing photon with 
the incoming lepton line in V\ ,2(0) we need to constrain 
y according to 



da 



, (2tt) 4 
~4fc ■ P l 
d 3 k' 



T^S^ik + P-k' -P' -q'y, 



d 3 P> 



d 3 q' 



2/^(2tt) 3 2P^(2tt) 3 2^(2tt) s 



(3) 



and after integrating out some of the phase space we are 
left with a five-fold differential cross section: 



da 



XbjV 



dx b] dyd\t\d(t>dip 16 7T 2 Q 2 v / ITe 2 



(4) 



The square of the amplitude receives contributions 
from pure DVCS (Fig. la), from pure BH (Figs, lb, 
lc) and from their interference (with a sign governed by 
the lepton charge), 



\r\ 2 = \T DVC sV+i+\T BH \ 

where the individual terms are given by 

„6 



(5) 
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(6) 



I^DVCsl'' 



DVCS 



y 2 Q 2 



1 = 



J2[c™ CS cos(n0) 

n=l 

±e 6 



„DVCS 



sin(n(/))l (7) 



x bj yHVi{(j))V2{4') 



(8) 



where the +/— sign in the interference stands for a neg- 
atively/positively charged lepton. 



V < Veal 



Q 2 + t 
Q 2 + x bj t 



(9) 



in order to avoid an artificially enhanced BH contribu- 
tion. This limit is only of practical relevance for fixed 
target experiments at very low energies. Collider exper- 
iments do not have any meaningful statistics for exclu- 
sive processes at very large y ~ 0(1). In the following 
discussion, I will neglect the contributions to the DVCS 
cross section containing transversity and terms higher 
than twist-3. 

The DVCS observables I will deal with later on are 
based on a less differential cross section than the five- 
fold one in Eq. (0J. The reason for this is first that the 
cross section in Eq. is frame dependent since the az- 
imuthal angles </> and tp are not Lorentz invariants and 
hence, they will be integrated out. Secondly, since a t- 
distribution is notoriously hard to measure, we also inte- 
grate out t, however with experimentally sensible cuts as 
will be discussed later. In consequence, our observables 
will be based on only a two-fold differential cross section. 
Note that the DVCS data currently available is at most 
for two-fold quantities, normally just one-fold or even 
totally integrated over. One might argue that the more 
variables in an observable are integrated out the more in- 
formation is lost, especially when studying higher twist 
effects. This is indeed true, however, one has to make 
a sensible compromise between wishful thinking on the 
one hand and experimental facts on the other. Also, as I 
will show below, these two-fold quantities are enough to 
clearly demonstrate the size of the higher twist effects on 
DVCS observables. In the following, I will concentrate 
both for the sake of brevity and the fact that these quan- 
tities are the easiest once to study, on the Single Spin 
Asymmetry (SSA) and the Charge Asymmetry (CA) de- 
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fined in accordance with experiments, the following way: 
2 f^d<p sin(</>)(A7 T -da 1 ) 



SSA = 



J Q 27r dcj) (da^ 



■dai) 



CA 



2 d(j) cos(4>)(da + — da ) 



(10) 



(11) 



/ c , # (da+ +da-) 

Here do^ and da^ refer to the two fold differential cross 
sections da/dxbjdQ 2 with the lepton polarized along or 
against its direction of motion, respectively; da + and 
da~ are the unpolarized differential cross sections for 
positrons and electrons, respectively. 

Even though I am trying to discuss the charge asym- 
metry (CA) for two experiments, EIC and CLAS, which 
cannot measure it at all since they are or will be running 
with electrons only, there exist experimental problems 
in measuring the proper quantity, the azimuthal angle 
asymmetry or AAA. The AAA is defined below 

r^/ 2 a„bh\ r 37r / 2 



AAA 



IZ% d^da - da BH ) - J:J' z d^jda - da BH ) 
Jq™ dtpda 



(12) 



where da BH refers only to the pure BH cross section. 
The experimental problem or challenge with the AAA is 
that it requires either a very good detector resolution i.e. 
many bins in <f> or an event by event reconstruction of the 
scattering planes. The last statement needs a word of ex- 
planation: Eq. l]12|l is equivalent to taking the difference 
between the number of DVCS minus BH events where the 
real 7 is above the electron scattering plane and where 
it is below that plane, divided by the total number of 
events. This procedure ensures that the numerator is not 
contaminated by BH, which would spoil an unambiguous 
interpretation of the observable in terms of the real part 
of DVCS amplitudes. Also, the only difference between 
Eq. Ijl 1(1 and 1(12(1 is due to the additional interference 
term in the denominator of Eq. 1(12(1 , which is small com- 
pared to the leading contribution, and a twist-2xtwist-3 
term in Eq. (JJJ), which is both suppressed in the kinemat- 
ics considered and small in the employed GPD model. 

Unfortunately, in the case of CLAS where BH is by 
far the dominant contribution in the cross section, one 
would need to subtract two large numbers for the above 
plane and below plane events inducing a huge statisti- 
cal uncertainty. Therefore, I will not discuss the CA for 
CLAS kinematics. 

Let me say a word, about the expected effects of 
uncalculated higher twist contributions besides the cal- 
culable twist-3 contributions. Assume an asymmetry 
A = xXy wnere B an( l X stand for leading twist 
contributions and and C and Y for the higher twist 
contributions in the interference and the cross section 
term respectively. Expanding the denominator yields 
A = B/X + C/X-B-Y/X 2 -C-Y/X 2 + . ... This shows 
that the leading higher twist contributions in A will orig- 
inate from the twist corrections to the leading twist in- 
terference term and will thus be of 0(y^-t/Q,rriN /Q) 



with respect to the leading term. Let us give a rough 
numerical example to see the significance of the higher 
twist corrections: Assume B = 1, X = 5,C = ±0.3 and 
Y = ±1, then A = 0.2 in the leading twist approxima- 
tion and A ~ 0.217 in the full result. This shows that 
even though we have 20 — 30% higher twist effects in the 
individual terms, they effectively cancel in the asymme- 
try if the corrections are both positive. However, if one 
or both of the higher twist corrections are negative, the 
result can vary between O(20 — 40%) from the leading 
twist result. The relative sign of the higher tiwst terms 
in both numerator and denominator will vary depending 
on the kinematic region one is exploring and, therefore, 
it is a priori not clear what the size of the higher twist 
corrections will be. Of course, these arguments only hold 
if the higher twist contributions do not become of order 
of the leading twist corrections. 

In Sec. [V] I will also talk about the one-photon cross 
section a(^*p) at small x^j defined through 



d 2 a(ep — > ep7) 
dydQ 2 

where 



= L a DV cs(l*P — > TP) 

ae.mjl + (I - y) 2 ) 
2nyQ 2 



with 



<tdvcs(i*P ->• IP) 



a 2 x 2 TT 2 

wes \t-0, 



Q 4 B 



(13) 



(14) 



and where B stems from the t-integration and will depend 
on both our cut-off in t and the model of the i-dependence 
I will choose for the GPDs. Furthermore, all higher twist 
effects are neglected in this quantity. 



III. THE GPD MODEL AND DVCS 
AMPLITUDES: TWIST-2 AND TWIST-3 

A. Modeling Twist-2 GPDs 

In the following I will use and review the model for 
twist-2 GPDs first introduced in [32|. 

Based on the aligned jet model (AJM) (see for example 
[H) the key Ansatz of ^ in the DGLAP region is: 



H 



S,V,g 



( x-C/A 

V 1-C/2 J 



l-C/2 



(15) 



where q l refers to any forward distribution and 
H S ' V [X,Q = H q {X,Q ± H^(X,Q. This Ansatz in 
the DGLAP region corresponds to a double distribution 
model I35L l3o | with an extremal profile function al- 
lowing no additional skewdness save for the kinemati- 
cal one. It will also be used for H and E. I will talk 
more about the exact details for H and E below. Note 
that I choose a GPD representation first introduced in 
0, which is maximally close to the inclusive case i.e 
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I e [0,1], ( = Xbj with the partonic or DGLAP re- 
gion in [Q, 1] and the distributional ampitude or ERBL 
region in [0, (]■ 

The prescription in Eq. (|15fl does not dictate what to 
do in the ERBL region, which does not have a forward 
analog. The GPDs have to be continuous through the 
point X — C and should have the correct symmetries 
around the midpoint of the ERBL region. They are also 
required to satisfy the requirements of polynomiality: 



M N = 



1 dXX 



N-l 



C 2-C 
(! + (-!)*% 



m(x,o-(-i) N - 1 H^(x,o 



XH«(X,Q 



+ (-1) 



A' 



c dXX N ~ 1 
2-C 



N/2 , . 2k 

I o 7 ) C 2k, \ 



k=0 



2-C 



W(X,C)+XHS(X,() 

(16) 



with X 



jzzff§ ■ The ERBL region is therefore modeled 
with these natural features in mind. One demands that 
the resultant GPDs reproduce the first moment Ma = 3 
and the second moment M 2 = 1 + CC 2 /( 2 ~ C) 2 JII- C 
was computed in the chiral-quark-soliton model |29| and 
found to be —3.2 and is related to the D-term [37j which 
lives exclusively in the ERBL region. This reasoning sug- 
gests the following simple analytical form for the ERBL 
region (X < Q: 

H g - v (X, C) = H a - V (C) [l + A 9 - v (C)C 9 - V {X, 0] , 

H S {X, C) = H S (C) ( ^fi) [1 + A S (C)C S (X, 0} , 



C/2 



(17) 



where the functions 



C* V {X,Q = - 



3 2 - C 



c s (x,c) 



2 C 
15 {2-C 



1 - 



X-C/2" 
C/2 , 



x_c/r 

C/2 , 



. (18) 



vanish at X = C to guarantee continuity of the GPDs. 
The A l (C) are then calculated for each C by demanding 
that the first two moments of the GPDs are explicitly sat- 
isfied. For the second moment, what is done in practice is 
to set the D-term to zero and demand that for each flavor 
the whole integral over the GPD is equal to the whole in- 
tegral over the forward input PDF without the shift. For 
the final GPD, of course, the D-term is added to the 
quark-singlet (there is no D-term in the non-singlet sec- 
tor) using the results from the chiral-quark-soliton model 
|29j. The gluonic D-term, about which nothing is known 
save its symmetry, is set to zero for Qq. Due to the 
gluon-quark mixing in the singlet channel, there will be a 
gluonic D-term generated through evolution. I will come 



back to this question in Sec. IIVI when discussing DVCS 
for the HERMES experimental setting. 

It would be straightforward to extend this algorithm 
to satisfy polynomiality to arbitrary accuracy by writing 
the A l (C) explicitly as a polynomial in £ where the first 
few coefficients are set by the first two moments and the 
other coefficients are then either determined by the ar- 
bitrary functional form, as is done here, or, perhaps the- 
oretically more appealing, one chooses orthogonal poly- 
nomials, such as Gegenbauer polynomials, for which one 
can set the unknown higher moments equal zero. Phe- 
nomenologically speaking, the difference between the two 
choices is negligible. 

The above Ansatz also satisfies the required positivity 
conditions |35l l38l l39l | and is in general extremely flexible 
both in its implementation and adaption to either other 
forward PDFs or other functional forms in the ERBL re- 
gion. Therefore it can be easily incorporated into a fitting 
procedure making it phenomenologically very useful. In 
what follows we will use MRST2001 and CTEQ6 
as the forward distributions for both LO and NLO. 

Let me quickly explain, why I only model certain 
C even and odd distributions in the quark sector, 
H S ' V (X, C) = Hi(X, C) ± H^(X, C). As one can see be- 
low, only the quark charge weighted H appears in the 
DVCS amplitude and not H v due to the C-even nature 
of the amplitude. The evolution equations for the GPDs 
are defined for the following C-even and C-odd singlet (s) 
and non-singlet (ns) flavor combinations: 



F i 

H v = H'l s = H q - - iJ ? , 
H s = — y^(H q + H q ) = —\^H s , 



N F 



(19) 



where H s mixes, of course, with the gluon and Hjf mix 
neither with each other nor with the singlet and the 
gluon. A single quark species, i.e. quark or anti-quark in 
the DGLAP region, or just a singlet or non-singlet quark 
combination in the ERBL region (due to the symmetry 
relations between quark and anti-quark in the ERBL re- 
gion in the off-diagonal representation of |j) can be ex- 
tracted the following way: 



i 



NS _|_ jjNS _ 



H s ) 



(20) 



in the DGLAP region and 

H s = {H^ s + H s ) , and H v = H n J , (21) 



in the ERBL region. Eqs. I|19I2UI21|) demonstrate that it 
is enough to model H S ' V in order to properly do evolution 
and extract the quark combinations relevant for DVCS. 

The construction of H proceeds analogous to that of 
H with opposite symmetries in the quark and gluon sec- 
tor and using the standard GRSV scenario |42| as the 
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forward input. Due to the change of symmetry in the 
ERBL region the analytical form changes to 

H S (X, C) = H S (C) [1 + A S (QC S {X, 0} , 
'X-t/2" 



H 9 ' V (X,C) = H 9 ' V {() 



[1 + A"> V {QC* V {X,Q] 



C/2 



(22) 



3 2 - C 



2 C 



C 9 ' V (X,() = 4 



2-C 

c 



1 - 



A -C/2 
C/2 

x-a/2 

C/2 




(23) 



Note that there is no D-term for the polarized GPDs due 
to symmetry requirements. 

For E, the asymptotic pion distribution amplitude i.e. 
the same Ansatz as in [23, was used (see also refer- 
ences therein for the same or similar Ansatze). 

The Ansatz for E used in this paper deviates from the 
one used in other studies to also include the gluon. First 
let me say a few general things about E: The symme- 
tries for E are the same as the ones for H in the ERBL 
region. Furthermore, one would naively expect that E as 
a function of C = Xbj dies out as C decreases i.e. behaves 
like a valence quark distribution. This is nothing but the 
statement that it becomes increasingly difficult to flip the 
spin of the proton as Xbj is decreasing. Also, we know 
that the first moment of E in the proton/neutron has to 
reproduce the respective anomalous magnetic moments 
(see for example This leads to 



k u — *2iKp ~\~ K n — 1.6T3 
t\>d — fcp ~\~ ^fcn — 2.033 . 



(24) 



Following [ij], the "forward" quark distributions from 
which to start is chosen to be 

1 



E u (x) = -u val (x) ■ K u 

E d (x) = d va i(x) ■ n d 
E s (x) = 



(25) 



For the DGLAP region and the quark singlet channel 
I will apply to Eq. (|25|) the same shift as in Eq. 1)150. 
The sea contribution is set to zero in the DGLAP region 
as was also done in ^3|. In contrast to |4i|, I choose 
to include a gluonic contribution. This contribution is 
important since these distributions are evolved from a 
low scale Qo to the relevant experimental scale in contrast 
to other groups [2^, which chose to neglect evolution 
effects. I model the gluon E 9 in the DGLAP region in 
the following way: 

First, we know that the total angular momentum of 
the proton J p for C = and t — is given by: 



J P = \=\ f dX [X(H S (X) + E S (X)) 



H 9 (X) + E 9 {X) 



(26) 



Since together with Eq. I|25|) I know now all the functions 
in Eq. I)2fi|l except the last, the contribution of the gluon 
to JP can be defined to be J 9 E = ±(1 - J s - J 9 H ). Thus 
one can define E 9 in the DGLAP region through 



E 9 (X,C) 



J 9 

E u G 

H 



(27) 



Note that in this particular model for the E q , the second 
moment of E s is zero due to symmetry, and since the 
second moment of the H 1 already saturates the angular 
momentum sum rule, the gluon contribution should be 
strictly zero. However, numerically, the second moment 
of the H l is never exactly 1, yielding a small but non-zero 
E 9 and secondly, the above construction is general and 
does not depend on the particulars of the model. 

In the ERBL region, I will use the same strategy as 
for H except that I require the sum rule for J p to be ful- 
filled after the shift in X rather than the first and second 
moment of just the H l . Keeping the symmetry require- 
ments for the E l in mind, this means that one can use 
Eq. (|221> and for the E i just with different A l {Q. 
In consequence, the shifted version of Eq. 125|) together 
with Eqs. Q22[l. I|23|) and J57J), gives a complete parame- 
terization of the E 1 . Note that this parameterization of 
the E l is in stark contrast to where the E l are of the 
same size or even larger than the H l even at small Xbj 
where this relation cannot hold. Hence the results for the 
asymmetries in this paper will vary from those in 28] at 
small Xbj whereas in the valence region they will be sim- 
ilar. Note furthermore that E contains a D-term in both 
the quark-singlet and the gluon. This D-term is identical 
to the one in H but enters with the opposite sign and 
thus cancels when considering the moments of the sum 
of H and E as done for the total angular momentum J p . 

As far as the t-dependence is concerned, I choose to 
model it the same way as in [2(| i.e. using a factorized 
Ansatz for the i-dependence from the rest of the GPD. I 
want to stress here that this is not really realistic theoret- 
ical assumption especially at larger values of t 01 ■ O n 
the experimental side, it was shown in |32| that in order 
to describe the ZEUS data on DVCS gH at large Q 2 , 
a Q 2 dependent slope of the f-dependence was required 
to describe the data. This in turn implies that the basic 
assumption of a factorized t dependence as well as the 
assumption that the t dependence of quarks and gluons 
is the same is wrong, at least at large Q 2 , since the Q 2 
dependence of the GPD can only be generated through 
perturbative evolution. The HI data on DVCS 0], for 
example, which lies in a lower Q 2 range does not a priori 
require a Q 2 dependent slope. This in turn means that 
at low Q 2 and low t, where most of the experimental 
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data lies, a factorized t dependence can still be used at 
the moment. The situation improves even more when one 
considers asymmetries since there either the t dependence 
partially cancels between numerator and denominator, if 
DVCS is dominant, or, if BH is dominant, the Q 2 range 
is such that a factorized approach is still not totally un- 
reasonable. As the accuracy and the kinematic reach of 
the data improves, however, one has to seriously address 
the issue of a non- factorized ^dependence of the GPD. 
I will discuss the issue in detail in Sec. Ivland propose a 
phcnomcnological solution. 

After having evolved the GPDs using the same pro- 
gram successfully employed in [3(|, the real and imagi- 
nary part of the twist-2 DVCS amplitude in LO and NLO 
give n below are calculated using the same program as in 

u 



«5(C,m 2 ,Q 2 ) = £< 
i 



2-C 

c 



2X 



dx T S a ,V/A l?±_ 1+i *\ jrS a ,V/A {x ^ c 3) 

V C M 
' 2X_ 

TV. 

2 



F J dX T S «- V ' A { I 



T g ' V/A (C a 1 O z ) - — 
1 dvcs^> r > I ~ at 



1 /2-C 



DVCS 
1 



N 



c 



dX T* V ' A f^-l + ie, %) ^ A (X, C, ,, 2 ) 
o V C M , 



± J dX T 9 > V I A | I 



2A Q 2 
TV 



(28) 



V/A stands for the vector/axial-vector i.e. unpolar- 
ized/polarized case and T stands for the appropriate 
GPD H, H, E or E. 

The +ie prescription is implemented using the Cauchy 
principal value prescription ("P.V.") through the follow- 
ing algorithm: 

rl r 2X 



C 



P.V.J dX T^—-1)F(X,(,Q Z 

rC /OX 

dX T [ — 1 ) (T(X, C, Q 2 ) - H(, C, Q 2 )) 



o 



J' dXT^- lj (T(X, C, Q 2 ) - C, Q 2 )) 



+ / dXT 



2X 

T 



- 1 



(29) 



The relevant LO and NLO coefficient functions can be 
found in Elllil. 



B. Modeling Twist-3 GPDs 

After having modeled the twist-2 sector which auto- 
matically takes care of the kinematic power corrections in 



the DVCS cross section as well, only the genuine twist-3 
sector remains. As shown in |20ll2lLl22| the twist-3 GPDs 
and thus twist-3 DVCS amplitudes can be expressed in 
the WW approximation through a combination of twist-2 
GPDs convoluted with a twist-3 coefficient function. 

The general structure of the twist-3 DVCS amplitudes 
can be found in Eq. (84) of [28( and reads in the repre- 
sentation of 



T tw - s (x b j,Q 2 )=T tw - 2 (x bj ,Q 2 ) 

4- n-u ■ 9 r V A (9) T ' 2m p x bj 

+ T^(x b ,,Q 2 ) 1 



( 1 %bj ) \ t ^mm ) 



T X (xbj,Q 2 ) 
(30) 



where the first three terms are the WW terms and the 
last term is a genuinely new dynamical contribution aris- 
ing from qgq correlations which will be neglected in the 
following. The T 1 - (x bj ,Q 2 )'s are linear combinations of 
the type C t ^,_ 3 <g> T which can be found in Eq. (87) of 
pEjj and do not have to be repeated here. Note that 
Eq. (H3J differs from Eq. (84) in 2§| by a factor (2 - x bj ) 
which I have pulled out for convenience. The convolu- 
tion C^_ 3 ® J 7 is done using Eqs. f2"%|) and f2l?)) with the 



C^_ 3 in LO having the following form 



2X 



In (1 - X/C + ie) 



(31) 



The respective subtraction factors P i.e. 
Jq dX T(^-—lj, for the imaginary and real part 
of the amplitude read 



fi e /^(0 = -C--CLi 2 (0- 7 Cln(C)^ 



Im 7^(0=--Cln(C). 



(32) 



q-tw-3 j e m particular the convolution and the deriva- 
tive of the convolution with respect to x b j , was computed 
numerically using an extended version of the program 
from 31] which will soon be available at 47]. 

In the WW approximation, one can use evolution of 
the twist-2 GPDs to evaluate Eq. H3(J|I at a scale different 
than the initial scale. This will allow one to study twist-3 
effects for the first time with a varying scale Q 2 . Since the 
twist-3 NLO coefficient functions are unknown (see [24[ 
for a recent calculation of the non-singlet quark sector) , I 
will restrict myself to a LO analysis in the twist-3 sector. 

Since twist-3 is entirely expressible through the twist-2 
in the WW approximation, I will use the same t depen- 
dence for the twist-3 DVCS amplitudes as I used in the 
twist-2 case. Having completed the specifications of the 
twist-2 and twist-3 sectors, I can now move on and dis- 
cuss twist-3 effects in DVCS observables which will be 
done in the next section. 
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IV. DVCS OBSERVABLES: TWIST-2 VS. 
TWIST-3 RESULTS 

In the following, I will discuss twist-3 effects in the SSA 
for four experimental settings: HERA, EIC, HERMES 
and CLAS and the CA for HERA, EIC and HERMES. 
Note that for illustrative purposes I will not only discuss 
twist-3 effects in LO but include the LO twist-3 ampli- 
tudes together with the NLO twist-2 and NLO kinematic 
power corrections terms. This is actually not legitimate 
since one is mixing different orders of a s . However, it 
serves both to illustrate the LO vs. NLO effects without 
genuine twist-3 effects and to set an upper limit on the 
twist-3 corrections in NLO, since the twist-3 NLO coef- 
ficient functions will not induce larger corrections than 
NLO in the twist-2 sector. This can be seen from the 
LO twist-3 coefficient function Eq. (|31() which has only 
a regulated logarithmic singularity instead of a regulated 
simple pole as the LO twist-2 coefficient function (see for 
example [2^, |^|). Hence it can be expected that the 
singularity structure of the NLO twist-3 will also be less 
severe than in the twist-2 coefficient function and the 
LO twist-3 effects will give a reliable upper bound for 
the NLO case. 



A. HERA 

In this section, I discuss the effects of LO twist-3 ef- 
fects on the CA and SSA in HERA kinematics with 
v 7 * = 319 GeV i.e. 27.6 GeV unpolarized/polarized 
positrons/electrons and 920 GeV unpolarized protons. 
Since it will be difficult for either ZEUS or HI to mea- 
sure a t distribution for DVCS, I will only discuss the 
CA and SSA integrated over t. Since the largest t for 
which the HERA experiments still have a signal is about 
~ — 1 GeV 2 , I choose a very conservative cut-off in t of 
—0.5 GeV 2 . I have checked that changing the cut-off to 
~ — 1 GeV 2 only alters the absolute answers on the order 
of 10% as well as leaving the relative twist-3 effect un- 
changed. Since I do neither know the acceptance curve 
in t for the HI and ZEUS detector, which induces an ad- 
ditional uncertainty in the answer, nor for any other of 
the experiments, for that matter, the chosen cut-off in t 
seems to be a sensible choice. 

As can directly be seen from the Figs. |3] and 0] the 
twist-3 effect in the CA are entirely negligible at HERA 
for both the MRST2001 and CTEQ6L LO parameteriza- 
tions. Furthermore, the two distributions give the same 
answer within about 10%. However, when comparing 
the respective NLO curves one finds differences of up to 
100%. They start to disappear for Xbj — > 0.1 in the given 
Q 2 range. For small Xbj, however, this difference disap- 
pears only for Q 2 > 40 GeV 2 . This feature was already 
noted earlier |26J in a pure twist-2 analysis and shown to 
be attributable to the very different NLO gluon distribu- 
tions at Qo. Also note that the exclusion of kinematic 
power corrections in ptjj lead to negative numbers for the 
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FIG. 3: t integrated CA in HERA kinematics vs. Xbj for two 
typical values of Q 2 and t max 
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FIG. 4: t integrated CA in HERA kinematics vs. Q for two 
typical values of Xbj and t max = —0.5 GeV 2 . 
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CA in HERA kinematics in contrast to our findings here. 
This illustrates the importance of kinematic power cor- 
rections for the CA. The NLO corrections, in particular 
at the smallest Xbj , are very large and only reduce for the 
largest Q 2 to about 100%. Again the same was found in 
a pure twist- 2 NLO analysis j2g and attributed to a large 
NLO gluon contribution in the real part of DVCS ampli- 
tudes. However, for Xb 3 > 10~ 3 and Q 2 > 10 GeV 2 the 
NLO corrections for CTEQ6M seem to be much smaller 
than in the case of MRST2001. I will come back to this 
point when I discuss the SSA. One word has to be said 
about the influence of the D-tcrm on the CA at this point. 
Based on the findings about the DVCS amplitudes in [3l| 
and by explicit comparison of results with and without 
a D-term, I conclude that the influence of the D-term is 
totally negligible for HERA. 

Turning now to the SSA, that we can see from Figs. [5] 
and that the twist-3 effects are even smaller than in 
the case of the CA. We also see that the room for twist-3 
effects in NLO is further reduced compared to the CA. 
Note that I discuss a positron rather than an electron 
beam and therefore the sign of the asymmetry is negative. 

Furthermore, we see that the NLO corrections are typ- 
ically of the order of 10 - 15% but at most 50%. This 
is in agreement with the results found in |26| | in a pure 
twist-2 analysis demonstrating that in contrast to the 
CA, the SSA is quite insensitive to kinematic power cor- 
rections. We also see that in NLO both MRST2001 and 
CTEQ6M give almost identical results, however, differ in 
LO for large Q 2 and Xbj < 10~ 3 which simply illustrates 
the fact that the LO gluon is larger for CTEQ6L than 
for MRST2001. In LO, this difference can only manifests 
itself after a longer evolution path since the DVCS am- 
plitude contains only quarks at leading order. In NLO, 
where the gluon enters directly on the amplitude level, 
differences in the gluon manifest themselves at lower Q 2 
in quantities very sensitive to the gluon contribution, like 
the CA due to its proportionality to the real part. This 
is well represented when comparing the NLO results for 
CTEQ6M and MRST2001 in the CA and the SSA at very 
small xi,j. 

The results for both LO and NLO suggest that both 
the CA and SSA should be easily measurable with fairly 
high precision at both the HI and ZEUS experiment! 



B. EIC 

In its current design the electron-ion-collider (EIC) will 
collide 1 — 10 GeV electrons from a linear accelerator with 
100 — 250 GeV unpolarized/polarized protons and unpo- 
larized ions of up to 100 GeV. Note that the projected 
luminosity for one year at the EIC will be larger than 
for the entire HERA run, enabling high precision studies 
of DVCS. For the figures below I chose a center-of-mass 
energy of y/s = 63.25 GeV which corresponds to a 5 GeV 
electron beam and a 200 GeV proton beam as a sort of 
average setting for the machine. The x^ range will thus 
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FIG. 5: t integrated SSA in HERA kinematics vs. Xbj for two 
typical values of Q 2 and t max = —0.5 GeV 2 . W stands for 
"with" and WO stands for "without". 
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be between roughly 10~ 3 — 10 _1 . Naturally, the higher 
y/s the closer the kinematics will be to HERA and, thus, 
also the results. Since we want to investigate the Xbj re- 
gion between HERA and HERMES with an overlap to 
both experiments, we do not want to go to the highest 
energies save for cross-checking HERA results. Let me 
start my discussion with the CA once again and then 
move on to the SSA. Note that I will apply the same t 
cuts as in the case of HERA to be able to compare the 
two settings. Though the proton target can be polarized, 
I will not discuss such an observable since it was shown 
in [2(| that such observables are essentially zero for a 
collider setting. 

As can be seen from Figs. and |S1 the twist-3 effects 
in LO for the CA are at most 10%, except at the small- 
est Xbj and lowest Q 2 = 2 GeV 2 where they can reach 
around 20% and have to be taken into account when try- 
ing to extract twist-2 GPDs from the data. In NLO the 
corrections seem larger but remember that this is only 
an upper estimate of the actual twist-3 effects in NLO 
and thus they are not more than 35% at the lowest Xbj 
and Q 2 . It is more likely, however, that they will be of 
the same size as the LO result or even smaller. Note also 
that the twist-3 effects quickly vanish for larger x^ within 
the entire Q 2 interval, which is mainly a kinematical ef- 
fect rather than a dynamical one. The two distributions, 
CTEQ6 and MRST2001, give very similar numbers in 
LO and at NLO for larger Q 2 , however, differ strongly 
at low Q 2 as already seen for the HERA setting. The 
relative NLO corrections are again large and follow the 
same pattern for both sets as at HERA. The influence of 
the D-term on the CA is again negligible. 

When looking at the SSA in Figs. I§1 and llOl one notices 
that the twist-3 effects are basically zero as in the case of 
HERA and that the NLO corrections are very moderate 
and of the same size for both sets as in the HERA case. 
Hence they can be safely neglected in a GPD extraction. 
Note that the shape of the SSA in Xbj and Q 2 is the 
mirror of the one at HERA since the EIC uses an electron 
rather than a positron beam. 

In conclusion one can say that except at low Q 2 and 
the smallest Xbj in the CA or a similar asymmetry, the 
twist-3 effects can be safely neglected and that the size 
is basically the same as in the case of HERA. These are 
very encouraging signs that, together with the high lu- 
minosity, DVCS will be measured with high precision at 
the EIC. Therefore, we will be able to reliably extract 
the leading twist-2 GPD H with high precision in a very 
broad range of Xbj and Q 2 from the EIC DVCS data! 



C. HERMES 

In the following I will discuss the fixed-target exper- 
iment HERMES with a center-of-mass energy of *Js = 
7.2 GeV. This allows, broadly speaking, to access a re- 
gion in x^ of about 0.05 — 0.3 with Q 2 from 1 — 9 GeV 2 . 
HERMES uses the electron/positron beam from HERA 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 



CA,L0wTw-3 
CA,LOwoTw-3 
CA,NL0wTw-3 
CA,NLOwoTw-3 

MRST2001 
Q 2 = 4GeV 2 




CA,L0wTw-3 

CA.LOwoTw-3 
- CA,NL0wTw-3 
CA,NLOwoTw-3 

MRST2001 Q 2 = 10GeV 2 




0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 





CA,L0wTw-3 




PA 1 fl uin Tuu *i 

On.LU WU 1 W-0 


- 


CA,NL0wTw-3 " 




CA,NLOwoTw-3 _ 




CTEQ6 




V\ Q 2 = 4GeV 2 


r's 

















10"' 



10" 



CA,L0wTw-3 

CA.LOwoTw-3 
- CA,NL0wTw-3 
CA,NLOwoTw-3 

CTEQ6 Q 2 = 10GeV 2 




FIG. 7: t integrated CA in EIC kinematics vs. Xbj for two 
typical values of Q 2 and t max = —0.5 GeV 2 . W stands for 
"with" and WO stands for "without". 
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FIG. 9: t integrated SSA in EIC kinematics vs. Xbj for two 
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FIG. 10: t integrated SSA in EIC kinematics vs. Q for two 
typical values of Xbj and t max = —0.5 GeV 2 . 



with an energy of E = 27.6 GeV. The gas target can be 
either unpolarized/polarized protons or unpolarized nu- 
clei. I will not discuss observables with a polarized tar- 
get for HERMES, since there is no clear leading DVCS 
amplitude, such as H in the case of the CA and SSA, 
and thus the disentangling of the various contributing 
GPDs is supremely difficult. In order to allow a compar- 
ison with the collider experiments, I once more choose a 
tmax of —0.5 GeV 2 which is also not an unrealistic choice 
given the fact that the average t for HERMES is about 
-0.25 GeV 2 . 

The CA, as can be seen from Figs. Illland ll2l receives 
larger twist-3 corrections in LO than the CA at HERA 
or EIC. However, except, at the lowest values of Q 2 and 
smallest Xbj where they can be as large as factor of 4, the 
corrections are generally speaking 15% or less. Note that 
as Q 2 increases the twist-3 corrections rapidly disappear 
in both LO and NLO. The LO results between the two 
sets agree very nicely but there is quite a difference in 
NLO. The NLO corrections themselves are again quite 
large but not larger than at HERA or EIC. In fact, for 
larger Q 2 and larger x j they are quite small. Note that 
when averaging the LO and NLO results with kinematic 
power corrections for the CA for both sets over Q 2 and 
x j one obtains the same numbers as in |32| while the 
number for LO with full twist-3 is about 0.1 neatly in- 
terpolating between the LO and NLO result of 0.12 and 
0.09 respectively. This compares very favorably with the 
experimental HERMES result of 0.11 ± 0.04 ± 0.04 jH 
for (x) = 0.12, (Q 2 ) = 2.8 GeV 2 , (t) = -0.27 GeV 2 . One 
can see that the averaging process washes out any differ- 
ences between the two GPD sets, which, however, were 
not that tremendous to begin with. A word about the 
D-term and its influence on the CA is in order at this 
point. Recently the first lattice results on the 

coefficient of the D-term were obtained and found to dif- 
fer from the prediction of the chiral-quark-soliton model 
p9| quite substantially. The respective calculations were 
done at different normalization points (/i = 2 GeV for the 
lattice i.e. within HERMES kinematics, and /j, = 0.6 GeV 
for the chiral-quark-soliton model). Evolution itself can- 
not account for the observed difference of about a factor 
of 4. When studying the LO and NLO evolution of the 
D-term using the GRV98 PDF with the above GPD 
Ansatz one finds that in both LO and NLO the quark D- 
term is reduced by about 30% from the respective input 
scale of Q = 0.51 GeV (LO) and Q = 0.63 GeV (NLO) 
to Q = 2 GeV (the difference between LO and NLO is 
about 2%), leaving still a factor of about 3 between the 
two results modulo the uncertainty associated with a glu- 
onic D-term which migh, given the right sign and size, be 
able to account for the observed difference. GRV98 was 
used since the input scales are very close to the one of 
the chiral-quark-soliton model. When studying the im- 
portance of the D-term for DVCS observables I find that 
if the D-term were either omitted or its size reduced by a 
factor of about 3 — 4, the CA would become so small that 
it would no longer be in good agreement with the data. 
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However, using a different Ansatz for twist- 2 GPDs based 
on a double distribution model (see for example |3|) and 
neglecting evolution effects, the authors of |2a describe 
the CA without a D-term. It was shown in [32j, though, 
that this type of double distribution Ansatz as chosen in 
[28l | , cannot describe the DVCS data in either LO or NLO 
when evolution effects are taken into account! The sit- 
uation will unfortunately remain unresolved until better 
fixed target data will become available. 

When turning to the SSA in Figs. 1131 and 1141 one can 
see that the twist-3 effects in LO are at most 10% and 
that the NLO corrections are, as in the case of HERA 
and EIC, very moderate and at most about 35%. The 
results of the two sets in LO are virtually identical and 
still within 20% at NLO. When averaged over Q 2 and Xbj 
the results of the two sets do not differ any longer and 
reproduce the LO and NLO results of (22J —0.28 and 
—0.23 as they should since the model is the same. This 
again compares favorably with the experimental result 
of -0.21 ± 0.04 ± 0.04 52] for virtually the same average 
kinematics as the CA. 

In conclusion, one can say that higher twist effects can 
be neglected for the SSA at HERMES and thus it can 
serve as a tool for GPD extraction. The CA is much 
more sensitive to twist-3 effects, however, they are still 
small enough that they can be neglected given the accu- 
racy of the data, except for the lowest Q 2 and Xbj values. 
This implies that for about Q 2 > 2-2.5 GeV 2 the CA 
can also be used for a GPD extraction or at the very least 
as a cross check to fits from smaller Xbj and the HER- 
MES SSA. The GPD model used in this study already 
produces very favorable agreement with the SSA and CA 
data without resorting to a fit and can thus serve as a 
basis for a successful parameterization. 



D. CLAS 

The CLAS experiment is a fixed target experiment 
with very high luminosity but a low center of mass energy. 
I will first investigate an electron beam of E = 4.3 GeV 
and then one with E — 5.75 GeV corresponding to the 
energies at the first and second CLAS run respectively. 
Here I will concentrate on the SSA and omit the CA 
or a similar asymmetry due to the mentioned difficul- 
ties CLAS has or will have with these type of asymme- 
tries as explained in Sec. |Hj Also, I will only discuss 
the set of GPDs generated from MRST2001 since the 
Qo = 1 GeV is low enough, compared to the one from 
CTEQ6 of Qq = 1.3 GeV, to have a meaningful range in 
Xf,j and Q 2 for both CLAS settings. 

Let us start with the lower energy setting, where I have 
chosen a t max = —0.25 GeV to get as large a range in 
Xbj and Q 2 as possible. As one can see from Fig. [TBI the 
twist-3 effects in LO are even for such a low energy as 
CLAS has less than 10% and thus basically negligible! 
This fits in nicely with the measured twist-3 effect at 
CLAS H which is about 10% of the measured SSA for 
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FIG. 11: t integrated CA in HERMES kinematics vs. Xbj for 
two typical values of Q 2 and t max = —0.5 GeV 2 . W stands 
for "with" and WO stands for "without". 
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FIG. 12: t integrated CA in HERMES kinematics vs. Q 2 for 
two typical values of Xbj and t max = —0.5 GeV 2 . 
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FIG. 13: t integrated SSA in HERMES kinematics vs. x bj for 
two typical values of Q 2 and t max = —0.5 GeV 2 . W stands 
for "with" and WO stands for "without". 



FIG. 15: t integrated SSA in CLAS kinematics vs. Xbj for 
four typical values of Q 2 and t max = —0.25 GeV 2 . W stands 
for "with" and WO stands for "without". 
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FIG. 14: t integrated SSA in HERMES kinematics vs. Q 2 for 
two typical values of Xbj and t m ax = — 0.5 GeV 2 . 



the central value but is compatible with zero within the 
experimental errors. Furthermore, the NLO effects are at 
most 50% and typically around 20% and thus not as large 
as one might have feared for such low Q 2 values. This is 
mainly due to the fact that the influence of the gluon on 
the amplitude in NLO at large Xbj is not as pronounced 
as at smaller x^j where its importance grows quickly. 
Notwithstanding this fact, the usage of perturbation the- 
ory at such small Q 2 remains still questionable on general 
grounds. However, one can definitely say that the twist-2 
handbag contribution to DVCS is the leading contribu- 
tion to the SSA at CLAS. This alone is quite an amazing 
statement given that one would naively have expected 
that at these energies higher twist contributions would be 
the dominant ones. When averaging over Xbj and Q 2 one 
obtains a value for the SSA in average CLAS kinemat- 
ics ((x) = 0.19, (Q 2 ) = 1.31 GeV 2 , (t) = -0.19 GeV 2 ) of 
about 0.2 in LO and about 0.14 in NLO which is, at least 
in LO, in good agreement with the experimental value of 
0.202 ±0.021 ±0.02 0. 

In the higher energy setting, I have introduced two dif- 
ferent t max values, in order to both compare to the lower 
energy setting and demonstrate how the SSA changes 
for a drastic change in cut-off for t. As can be seen from 
Figs. 1161 1171 118l andll9lthe twist-3 contributions are again 
very small in both LO and NLO for both cuts in t and 
are at most 10% which is in agreement with the value 
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at lower energy. The NLO corrections are mostly mod- 
erate except for the lowest values of x^j and Q 2 as seen 
at lower energies. The distribution in Xbj for different 
Q 2 between the lower and higher energy setting at the 
same t max (Figs. El and 1161) shows the distributions to 
be very similar both in shape and size. When compar- 
ing the different cuts in t for the higher energy setting 
(Figs. El an d El an d Figs. El an d I19|) one notices that 
the distributions in xi,j with the higher t cut are wider 
and thus flatter than the one for a lower cut in t. The 
maxima of the curves move also to larger values of Xbj- 
The Q 2 distributions are virtually unaltered. There is an 
overall tendency for the maxima to be somewhat higher 
for the higher t cut, but only by at most 10%. 

We thus see that different cuts in t have only a marginal 
effect in the size and distribution of the SSA. When fur- 
ther averaging over Xbj and Q 2 the sensitivity will be 
further reduced. In fact, as one can see from the figures, 
cuts in Xbj and Q 2 will have a much bigger effects than 
the one in tl 

Unfortunately, there is no published data yet with 
which to compare and without knowing the average kine- 
matics, let alone the experimental acceptance, it is im- 
possible to make a sensible prediction at this point. 



V. PHENOMENOLOGICAL 
PARAMETERIZATION OF THE GPD H 

Let me now turn to a sensible, phenomenological pa- 
rameterization of the X, ( and t dependence of the lead- 
ing twist-2 GPD H at a low normalization point Qq. 
First, I will give the parameterization of H in LO and 
NLO and then justify it based on the available data. 

As is clear from the preceding section, the NLO pa- 
ramcterizations seem to work very well in their current 
form except for CLAS. However, for CLAS the GPD H 
will not be the leading one anymore as it is for HERA, 
EIC and HERMES. In fact, the contributions from other 
GPDs could be set to zero without changing the HERA 
and only by a few percent the HERMES results. Since I 
only want to make a statement about H, I will restrict 
myself to a good description of the data from HI, ZEUS 
and HERMES. 

Based on the analysis carried out in 32], the 
MRST2001 NLO PDF parameterization at Q = 1 GeV 

with h-q'cD = 323 MeV 

using the prescriptions of 
Eqs. i|15|) . l|17fl an d 1|18|> does the best job in describing 
the DVCS data from HI, ZEUS and HERMES. There is 
no need to change the NLO parameterization of Sec. lIIII 
The story is different for LO. The LO results using 
Eqs. ifUfl) . (|T7jl and i(Tg|) are consistently above the DVCS 
data save for CLAS. A way to find a LO parameteriza- 
tion giving a good description of the data is to vary the 
shift parameter a in £. The shift parameter is given by 
the number in front of £ in the argument of the forward 
PDF (X - a()/(l - <) i.e. a = 1/2 in Sec.|Hl] A shift 
parameter of a = 1/2 works well for NLO but not for 
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LO. In fact, the best description of the DVCS data in 
LO is found for a ~ for the MRST2001 LO PDF with 
Qo = 1 GeV and Aq£,~ 4 ' LO = 220 MeV. Note that any 
a ^ 1/2 will violate the "Munich symmetry" of double 
distributions [35|, |53| though still fullfilling all the other 
requirements. Given the fact that the NLO parameter- 
ization fullfills all necessary requirements and, save for 
absolute numbers, looks the same as the LO parameter- 
ization, one can conclude that the LO parameterization 
is still phenomenologically useful since it describes the 
data and can be used for good quantitative estimates 
though it neglects higher order corrections and violates 
some subtle symmetries. 

Let me illustrate this with the example of the HI data 
Figs. 1^01 and [2] for the DVCS 7*-proton cross section, 
cr(7*p), Eq. lfH|) . As I explained before, the leading 
DVCS amplitude at small Xbj is generated via the GPD 
H . The interference term in the DVCS cross section is, 
after integration over 0, only a percent contribution to 
a(j*p). The BH term is usually also negligible compared 
to the pure DVCS term, however, since we are able to 
compute it unambiguously to high accuracy, one can sim- 
ply subtract the BH contribution from the data. In this 
case the t dependence can be simplified to an exponential 
form e Bt . For the HI data it is sufficient to take B to be 
a constant. However, for the ZEUS data this is not suf- 
ficient anymore (see [12 )• Though not necessary, I will 
use the Q 2 dependent slope of |32l | 



B(Q 2 )=B 1-Cln 
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Ql 



(33) 



with B = 8 GeV" 2 , Q = 2 GeV 2 , C = 0.15. The 
reason for choosing such a parameterization are given in 
great detail in [33] and need not be repeated here. A 
physically intuitive explanation for this behavior of the 
slope is given in 

As can be easily seen when comparing the upper and 
lower plots in Figs. H3 and E2 the LO MRST2001 curve 
now compares very favorably with the HI data. In fact, 
it gives virtually the same result as CTEQ6M i.e. it 
underestimates the ZEUS data somewhat. For the fixed 
target kinematics it is in agreement with the HERMES 
data on the SSA (0.21) and CA (0.09) and the CLAS 
data (0.17) on the SSA when kinematically averaged. 

I would like to comment now on the t dependence of 
GPD H at Qo- Note that in the parameterization of 
the slope B, I did not introduce a Xbj or W dependence 
as is customarily done (see for example and refer- 
ences therein) to account for cone shrinkage i.e. the fact 
that the slope increases as Xbj decreases for constant Q~ . 
However the slope change in x^j for HERA kinematics is 
only of the order of 10% and can thus be neglected for 
practical purposes. Furthermore, the necessity of a Q 2 
dependent slope signals a breakdown of factorizing the t 
from the X and £ dependence as has always been done 
in modeling GPDs. The breakdown at small Xbj does 
not occur until fairly large values of Q 2 which is very 
suggestive of the following scenario: At the initial scale 



Qo one has a factorized component of the t dependence 
which serves as a normalization and will be different for 
valence-quarks, sea-quarks and gluons. This difference 
in normalization between quarks and gluons will change 
as the gluon mixes with the quark-singlet under pertur- 
bative evolution. This change in the normalization will 
be Q 2 dependent and thus, the form of Eq. i|33|) is very 
natural since evolution resums logs of Q 2 . The Xbj de- 
pendence of the slope can be generated similarly if the £ 
dependence of the GPD is £~ A - Qt i. e . a regge-like depen- 
dence with a ~ 0.25 GeV -2 but possibly smaller, espe- 
cially for gluons, which will also acquire a logarithmic Q 2 
dependence through evolution much like the logarithmic 
slope of F2 54j . This type of slope change can also be 
parameterized with an exponential form as in Eq. I|33|l 
with Q 2 replaced by W i.e. Xbj- Note that in order to 
maintain polynomiality of the GPD H for t ^ 0, the coef- 
ficients A 1 (C) in Eq. (|17fl will also acquire a t dependence 
in order to compensate the extra factor of C,~ at . 

A sensible parameterization for the t dependence of 
the GPD H at Qo would thus be to choose a factorized 
exponential part with a square root of the slope of \B q ~ 
4-5 GeV" 1 for the quark sea and \B g ~ 2 - 2.5 GeV -1 
for the gluon. The valence quarks retain the dipole dis- 
tribution in t used in this paper. The t dependence in £ 
can then be chosen to be (~ at for the quark sea and the 
gluon with a ~ 0.25 GeV~ 2 for simplicity. The valence 
distribution could in principle also have a £ _Qt behavior 
but this will be extremely difficult to disentangle from 
data from small t. A more accurate parameterization in 
t will require much more precise data as can be expected 
from the EIC. 

This completes the phenomenological parameteriza- 
tion of the input GPD H in LO and NLO. 



VI. CONCLUSIONS 

I have given a detailed account of LO twist-3 effects 
in the WW approximation including their perturbative 
evolution on DVCS observables for kinematical settings 
equivalent to the HERA, EIC, HERMES and CLAS ex- 
periments. Based on the successful GPD Ansatz of [32$ . 
I found that the twist-3 effects for the collider settings 
are negligible save for the lowest values of Q 2 and Xbj- 
For these Q 2 and Xbj values the twist-3 effects still only 
reach about 10% in observables sensitive to the real part 
of DVCS amplitudes namely the charge asymmetry and 
even less in observables sensitive to the imaginary part 
of DVCS amplitudes such as the single spin asymmetry. 
The twist-3 effects for the fixed target experiments were 
only sizeable for the charge asymmetry at low Q 2 and 
Xbj, however not larger than 10— 15% for the single spin 
asymmetry. The common feature, of course, is the virtual 
disappearance of these effects for Q 2 values larger than 
about 3 — 5 GeV 2 depending on the value of x^. The 
relative smallness of twist-3 effects combined with the 
fact that twist-3 DVCS amplitudes in the WW approx- 
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imation are entirely expressible through twist-2 GPDs 
makes an extraction of at least the unpolarized twist-2 
GPD H which is leading in at least three of the four kine- 
matical settings, entirely feasible even with the current, 
relatively low statistics data. The EIC with its high lu- 
minosity will then enable a high precision extraction of 
the twist-2 GPD H. 

Since the current data from HI, ZEUS and HERMES 
are already remarkably restrictive for H, I give a first, 
phenomcnological, parameterization, though not a fit, in 
A, ( and i at a low normalization point Q which de- 
scribes all available DVCS data from HERA and HER- 
MES in both LO and NLO and from CLAS in LO. 
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